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Within the framework of Skyrme energy density formalism, we investigate the role of surface 
corrections on the fusion of colliding nuclei. For this, the coefficient of surface correction was 
varied between 1/36 and 4/36, and its impact was studied on about 180 reactions. Our detailed 
■ investigations indicate a linear relationship between the fusion barrier heights and strength of the 

surface corrections. Our analysis of the fusion barriers advocate the strength of surface correction 
^ ! of 1/36. 

Q\ • PACS numbers: 25.70.Jj, 24.10.-i. 



I. INTRODUCTION 

b _ 

The usefulness of the Skyrme energy density formalism (SEDF) in understanding the fusion dynamics has been 
i " i ' widely accepted 3[ . Several ion- ion potential models and parametrizations have been suggested based on the 
same formalism [lH3|. Interestingly, along with SEDF, many other models/formalisms based on either micro- 
scopic/macroscopic or phenomcnological picture are found to reproduce the fusion barrier heights within ±10%. 
In a recent study we employed as many as 16 versions of different proximity potentials and found that one can 
"^T , reproduce the barrier heights across the periodic table within ±8%. Even in the last few years, many sets of the 
^y-j ° Skyrme forces that reproduce the ground states properties of large number of nuclei have also been proposed [H-Q- 
^-j. . Within the same SEDF, the use of different Skyrme forces can yield difference of ±10% for fusion barriers. The same 
Skyrme forces have also been used at intermediate energies to investigate many rare phenomena [5(. The strength 
and form of nuclear potential is also very important in the study of cluster decay • 

Generally, SEDF consists of Hamiltonian, that depends on the nucleonic density, kinetic energy density as well as 
H ! on the spin density [l|. Among all these, the form and strength of the kinetic energy density has always controversial 
in the literature pj. Mostly reported works used approximations based on the Thomas-Fermi approximation pzl— Tiolj . 
. | In the literature, additional correction in term of gradient term has been suggested over and above the Thomas- 
Fermi approximation @-Qj]]. Though all studies advocate the inclusion of this term, its strength has yet not been 
resolved and many different strengths are available in the literature 0, EH EH . We plan to explore the role of this 
surface correction term in heavy-ion collisions via fusion process and want to understand whether one can narrow 
down the choice of this parameter or not. We shall also attempt to present a simple parametrization of this term for 
fusion barriers. This aim is achieved by employing SEDF within the proximity concept discussed in brief in section 
HTl The results are presented in section IIIII and a brief summary is presented in section IIVI 



II. SKYRME ENERGY DENSITY FORMALISM: 



In the Skyrme energy density formalism [l[, the nuclear part of the interaction potential Vn{t) is calculated as a 
difference of the energy expectation value at a distance r and at infinity (i.e. at r = oo): 

V N (r) =E{r)-E (oo) , (1) 



'Electronic address: idsharma.pu@gmail.com 



2 



with 

E = J H{r) dr. (2) 
In our formalism, the energy density functional H (r) read as; 

H(p,T,J) = A-r+it [(l + ix )p 2 -(x + i)(p2 + p 2 )] + l( tl +t 2 )pT 

+ o(*2 ~ k){PnT n + PpT p ) + —(t 2 - 3ii)pV 2 p 

8 16 
+^2"( 3 *i + i 2)(/O n V 2 p„ + p p V 2 p p ) + -hp n p p p 

--W (pV • J + p„V ■ J n + p p V • J p ). (3) 

Here p is the nucleon density taken to be two-parameter Fermi density and J is the spin density which was generalized 
by Puri et al., for spin- unsaturated nuclei The remaining term is the kinetic energy density r. The most widely 
used Skyrme force type SIII [l3[ is used in the present analysis. 

The evaluation of kinetic energy density term was done within the Thomas-Fermi (TF) approximation which is 
a well known alternative to the Hartree-Fock (HF) method. As shown by Myers and Swiatecki [8] and Hilf and 
Siissman [9], the kinetic energy density r can be separated into volume term To and surface term t\ plus reminder. 

— 5 

In the first order approximation, one can limit to t = I (§tt 2 ) 3 p 3 term only. However, due to the serious drawback 
in explaining the nuclear surface as well as the densities, a demand of further correction to this term was realized. 
Bethe and Brueckner (loj proposed that this difficulty can be overcome by adding some gradient term and correction 

proportional to ^—y—', derived long ago by von-Weizsacker [11]. Alternative to this is the gradient extension method 
that also yields similar results. The coefficient in the above term was first thought to be close to 1/4 ( = 9/36) [3], 
later on, it was found that this strength does not give proper results for kinetic energy density. Following this, several 
authors suggested different values for surface correction 

[Ml that ranges from 1/36 to 9/36. The kinetic energy 

density can therefore be written as; 

(vp) 2 

T = T + T A = T + A- (4) 

P 

where A is an adjustable parameter. As discussed by Graf 0, its value is supposed to be between ^ and Jg. In 
last few years, large number of theoretical models however used the value of A = 0, HH, El- One is wondering 
how A can alter the results of fusion barriers and further can we understand its effect in terms of some parameterized 
form or not. It is worth mentioning that different strengths of A can be thought to be close to Yukawa term used in 
addition to the Skyrme forces at intermediate energies ||. Therefore, this term is very important for reproducing the 
surface properties of nuclei as well as for multifragmentaion. We shall calculate the ion-ion potential using the above 
formalism for different values of A and then extract the fusion barrier heights and positions. 



III. RESULTS AND DISCUSSIONS 



In the present study, as many as 180 reactions involving even-even masses between 24 and 246 are taken for analysis. 
As noted, symmetric as well as asymmetric colliding nuclei have been taken into account. We calculated the ion- 
ion potential for all 180 reactions using different surface strengths. The heights and positions of the barriers were 
calculated using conditions: 

dV T {r) , d 2 V T (r) , 

— |r=fl B =0, and —— — r =R B < 0. (5) 

ar dr z 

For the present analysis, we took A = and The barrier height is denoted by Vb and its corresponding 

position is marked as Rb- We shall also attempt to parameterize the fusion barriers, thus, obtained and will present 
unique correlation between different values of A. 
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In Fig. 1, we display the percentage difference of the barrier heights and positions calculated using different values 
of A (between and ^) over experimental values defined as: 

ytheor _ yexpt 

^Vb % = - JL —^~ x 100, (6) 

V B 

as a function of Z\Z 2 - The experimental data are taken from Refs. 0, [l5 L 24]. We see that different values of A 
can explain barrier heights within ± 10% of the experimental values. Whereas, due to large uncertainty in the 
measurement of barrier positions the deviation in some cases can be quite large ± 30%. These data have been of 
controversy in the literature 1]. Due to different experimental setups, measurements do not yield few barrier positions 
as per known trend. A more careful look reveals that fusion barrier heights calculated with A = ^ yield better results 
on average as compared to other values of surface correction used. From this figure, which is spanned over the entire 
periodic table, it is evident that A = ^ may be preferred as the coefficient of the gradient term in the kinetic energy 
density correction. This finding is in agreement with other large number of studies based on the structural aspect of 
heavy-ion reactions 0, fl2l fl3 . However slight deviations in some cases may be due to the experimental uncertainties 
reported in many papers [241 ] or it may also be due to the influence of additional higher order correction terms that 
we have not taken into account in the present calculations. 

We shall now attempt to parametrize the fusion barrier heights and positions using different A - values in terms of 
charges and masses of the colliding pair. In Fig. 2, we display the barrier positions Rb as a function of {AX' 3 + A\ /3 ) 
and Vb as a function of z i^i (1~ p»»oi )i where R R nal is the analytical barrier positions obtained after parametrization 

B B 

have been performed over Rb results. When we increase the value of A from to and a monotonous 

increase in the nuclear part of the potential results in the reduction of the height of the barrier and therefore, the 
barrier positions are pushed outward. The addition of the strength A increases the attractive part of the interaction 
potential. In other words, one can counterbalance the repulsive Coulomb potential at larger distances, therefore, 
pushing the barrier outwards. As a result, net decrease in the barrier height occurs. In Fig. 2, we display the results 
with A = jjg, = Jg, and = The left part is for the barrier positions whereas right part is for the barrier heights. 
The barrier heights and positions are parametrized in terms of the following relations: 



R anai = a + b Xl] and V^ nal = c X 2 , (7) 

where X\ — (A^ 3 + A^ 3 ), and X 2 = Z ^nai (1 — R <Lai )■ Here a , b , and c are the constants displayed in the 
figure. We see a linear increase in the barrier heights and positions with the masses of the colliding nuclei. This is 
in agreement with many previous studies 25]. With the increase in the strength of A, a monotonous decrease in the 
barrier height and increase in the barrier positions can be seen. In all cases, a linear fit can explain the effect of A on 
the barrier positions very well. 

In Fig. 3, we display the percentage difference between the analytical and exact values as a function of Z\Z?,. Very 
interestingly, we note that we can reproduce the actual barriers (heights as well as positions) within ± 2.5% for all 
values of A. This introduces great simplification in calculating the barrier positions using different strengths A. It 
would be of further interest to understand how different values of A affects the fusion barriers in different reacting 
nuclei. 

In Fig. 4, we display the parametrized fusion barrier heights V B mal , as a function of the different values of ^ for 
the reactions of 24 Mg + 24 Mg, 28 Si + 28 Si, 48 Ca + 48 Ca, 64 Ni + 64 Ni, 40 Ca + 50 Ti, and 48 Ti + 58 Ni. Interestingly, 
we note that the barrier heights reduce systematically as one move from A = 0.0 to A = ^ values. This reduction 
can be further quantify by the relation: 

V «nal =a{1 _ § *_ ) (8) 

Here a and S are the constants whose values depend upon the colliding pair (see Fig. 4). Therefore, decrease in the 
barrier height with A is on expected lines, a linear reduction is of particular interest. 

Finally, in Fig. 5, we display the fusion cross sections <jf us calculated using the formula of Wong [26|, as a function 
of the center-of-mass energy, E c .rn, for the reactions of 24 Mg+ 24 Mg [MOlt 28 Si+ 28 Si [TM1, A8 Ca+ 48 Ca [UHJ 



64 Ni + 64 Ni [21|, 40 Ca + 50 Ti [23|, and 48 Ti + 58 Ni [22j. We see that no particular value of A explains the fusion 
cross sections. In few cases, however, A ~ ^ yields better results whereas in other cases, higher values of A has better 
edge. 
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IV. SUMMARY 



In the present study, we investigated the effect of surface corrections on the fusion process. Our finding 
over 180 reactions reveals that A = ^ can be a better choice for the surface correction and it is in agreement 
with earlier attempts. We also obtained parameterized form of the fusion barrier heights and positions for dif- 
ferent strength of A. We further found that the fusion barrier heights depend linearly on the strength of the coefficient. 
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FIG. 1: (Color online) The percentage difference between the theoretical and experimental values as a function of Z\Zi 
different values of surface correction A. 
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FIG. 2: The barrier positions Rb (left side) and barrier heights Vb (right side) as a function of Al^+Ay 3 and z ^nai (1~ naLai )i 
respectively for A = 1/36, A = 3/36, and A = 4/36. The solid lines represents the straight line linear fit over the points. 
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FIG. 3: The percentage difference between the analytical and theoretical values as a function of Z1Z2 for different values of 
surface correction A. 
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FIG. 5: (Color online) The fusion cross sections trj us (mb) as a function of the center-of-mass energy E c . m . for the reactions of 
24 Mg + 24 Mg, 28 Si + 28 Si, 48 Ca + 48 Ca, 64 Ni + 64 Ni, 40 Ca + 50 Ti, and 48 Ti + 58 Ni using different values of surface correction 
A. 



